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Abstract 

The effect of the " chromo-electric" dipole moment on the electric dipole mo- 
ment (EDM) of the neutron is studied in the two-Higgs-doublet model. The 
Weinberg's operator O^g = GGG and the operator Oqg = qaGq are both inves- 
tigated in the cases of tan/3 ^ 1, tan/3 <^ 1 and tan/3 ~ 1. The neutron EDM 
is considerably reduced due to the destructive contribution with two light Higgs 
scalars exchanges. 
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1 Introduction 



The electric dipole moment (EDM) of the neutron is of central importance to probe a 
new origin of CP violation, because it is very small in SM ~ 10"'^° — 10~^^e-cm). 

Begining with the papers of Weinberg there has been considerably renewed interest 
in the neutron EDM induced by CP violation of the neutral Higgs sector. Some studies 
P, ^ H revealed the importance of the " chromo-electric" dipole moment, which arises 
from the three-gluon operator GGG found by Weinberg [0 and the light quark operator 
qaCq introduced by Gunion and Wyler 0, in the neutral Higgs sector. Thus, it is 
important to study the effect of these operators systematically in the model beyond 
SM. We study the contribution of above two operators to the neutron EDM in the 
two-Higgs-doublet model(THDM) P]. The 3x3 mass matrix of the neutral Higgs 
scalars is carefully investigated in the typical three cases of tan (3 ^ 1, tan/? ~ 1 and 
tan j3 1. In this model CP symmetry is violated through the mixing among CP = + 
and CP = — Higgs scalar states. 

In order to give reliable predictions 0, one needs the improvement on the accuracy 
of the description of the strong-interaction hadronic effects. Chemtob P] proposed a 
systematic approach which gives the hadronic matrix elements of the higher-dimension 
operators involving the gluon fields. We employ his model to estimate the hadronic 
matrix elements of the operators. 

2 CP violation parameter in THDM 

The simplest extension of SM is the one with the two Higgs doublets |0. This model 
has the possibility of the soft CP violation in the neutral Higgs sector, which does not 
contribute to the flavor changing neutral current in the B, D and K meson decays. 
Weinberg has given the unitarity bounds for the dimensionless parameters of the 
CP nonconservation in THDM. However, the numerically estimated values of these 



parameters are not always close to the Weinberg's bounds 0. Although it is difficult 
to estimate the magnitudes of the CP violation parameters lmZi{i = 1,2) generally, 
we found that the neutral Higgs mass matrix is simplified in the extreme cases of 
tan/5 <^ 1, tan/3 ^ 1 and tan/3 ^ 1, in which the CP violation parameters are easily 
calculated. The CP violation parameters ImZ^^^ are deduced to 
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ImZ^ = — -tW2 u\ , 1 
sm/5 

where uf'^ denotes the i— th component of the k—ih. eigenvector of the 3x3 Higgs 
mass matrix and tan (3 = f2/t^i(t^i(2) is the vacuum expectation value of $5(2) giving 
the masses of (i(u)-quark sector). 

In this model, Higgs potential is generally given as 
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where the parameters satisfy the conditions |T0| 
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It is noted that, in the case of MSSM, SUSY imposes the conditions on the parameters 



9i = 92 = \{9w^ + 9w'^), 
9 = \{9w^ -9w'^), 



9 — —29w , 



h = 0. 
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Here h=0 means that in MSSM CP violation is not caused throygh Higgs sector. The 
simplest SUSY extention from MSSM that can have CP violation in the Higgs sector 



is also discussed HIT 



Let us estimate uf'^ by studying the Higgs mass matrix whose components are 



M- 



11 



33 



^12 



MI3 



2gi\vi\^ + 9'\v2\^ + 



{\V1\' + \V2\') 

1^1^21(2^7 + ^7') + 



^ + Re{hvfvl) 
^ + Rejhvfv^) 

9 + 1 ^ 

Re{hvfv^) - ^ 

\VlV2\ 



\VlV2 



lm{hvl'v'2) , 



+ \V2 



-lm{hvl'v'2) . 



As a phase convension, we take h to be real and 



(5) 



vl V2 = \vi\ \v2\ exp(2i0) . 



(6) 



At first, we consider the case of tan (3^1 with retaining the order of cos (3 and setting 
cos^/9 = and sin/3 = 1. Then, the mass matrix becomes simple, so the eigenvectors 



of in Eq.(5) are easily obtained as follows: 

= { cos /3 - e sin /?, -sin/?, 0}, (7) 
u*-^^ = {sin jSc^, (cos/5 — e sin/5)c0, — s</,}, 
u^^^ = {sin/5s0, (cos/? — e sin/3)S(/,, c^}, 

where c(s)<^ = cos(sin)0 and 

cos/7 . 8 

The diagonal masses are given as 

M^ = 2g,, Ml = g' + l+h, Ml = g' + 1 - h (9) 

in the v"^ = v'l + unit. The lightest Higgs scalar to yield CP violation is the second 
Higgs scalar with the mass M2 since ^ is positive from Eq.(3) and we take h to be 
negative as convention. The Higgs scalar with Mi does not contribute to CP violation 
because of v!-^^ = 0. The absolute values of g' is expected to be 0(1), but h seems to 



be small as estimated in some works [12, 0. For example Froggatt et al. give the 
numerical values for the parameters in the case of tan/? ^ 1 by using infrared fixed 
point analysis through the renormalization group equations as 

gi ~ 0.96, g2 ~ 0.88, g ~ 0.82 

g' ^ -1.20, h ~ -0.09. (10) 

Therefore, the masses M2 and M3 may be almost degenerated. Then, CP violation 
is reduced by the cancellation between the two different Higgs exchange contributions 
ImZj^^"* and ImZ^"* since ul^^l^^ and u^f\l^\i=l,2) have same magnitudes with op- 
posite signs. Thus, it is noted that the lightest single Higgs exchange approximation 
gives miss-leading of CP violation in the case of tan/? ^ 1. 

For ImZi, our result reaches the Weinberg bound, but for ImZ2 the our calculated 
value is suppressed compared with the Weinberg bound in the order of 1/ tan/3. 



CP violation in the case of tan/5 <^ 1 is similar to the one of tan/5 ^ 1. For 
lmZ2, our numerical result reaches the Weinberg bound, while for ImZi the calculated 
value is suppressed from the Weinberg bound in the order of tan /3. The relative sign 
between ImZi and lmZ2 is just the same as in the case of tan/5 ^ 1. 

The last case to be considered is of tan /? ~ 1. In this mass matrix, the off diagonal 
components are very small compared to the diagonal ones because gi ~ g2 is suggested 
by some analyses []T^ |T3] and h is also small as in the case of tan/? ^ 1. We can 



calculate ImZi by fixing both values of h and M2/M3. For both lmZ2 and ImZi, the 
calculating values are roughly 1/3 of the Weinberg bounds. The relative sign between 
ImZi and lmZ2 is opposite. 

3 Formulation of the neutron EDM 

The low energy CP- violating interaction is described by an effective Lagrangian, 

Lcp = Y.c^{M,^i)o,{^i) , (11) 

i 

where are the three gluon operator with the dimension six and the quark-gluon 
operator with the dimension five as follows: 

9s abc ~a b 

03g{x) = — "C^^^G^^Gl^ , (12) 

where q denotes u,d ot s quark. The QCD corrected coefficients Cj are given by the 
two- loop calculations in Refs. P, Bl. The coefficients Cj are given as 
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where the functions f{x),g{x),h{x) are deduced from loop integral as given in Refs. 



For the strong interaction hadronic effect, the systematic technique has been de- 
veloped by Chemtob |^ in the operator with the higher-dimension involving the gluon 
fields. The hadronic matrix elements of the two operators are approximated by the 
intermediate states with the single nucleon pole and the nucleon plus one pion. Then, 
the nucleon matrix elements are defined as 

{N{P)\Om\N{P)) = AdJ{P)ii,U{P), 
{N{P')\0,\N{P)Tx{k)) = BJJ{P')T-U{P) , (14) 

where U{P) is the normalized nucleon Dirac spinors with the four momuntum P. Using 
Ai and Bi{i = ug, dg, sg, 3g), the neutron EDM, d^, are written as 

= ^ E CiA. + F{9.nn) E C^B, , (15) 

where /i„ is the neutron anomalous magnetic moment. The -F( (7,7 at at) was given by 
calculating the pion and nucleon loop corrections using the chiral Lagrangian for Nt^'-j 
[^. The coefficients Ai and Bi were given by the large current algebra and the 7]q 
meson dominance P]. 

4 Numerical results of the neutron EDM 

Let us begin with discussing the numerical results in the case of tan/? ^ 1. The con- 
tributions of Oug and are are negligibly small because the CP violation parameters 
are roughly estimated as 

ImZf ^ ~ -ImZp^ ~ ^ tan^ p. (16) 



The main contribution follows from the one of Odg + Osg, in which the operator Ogg is 
dominant due to the s-quark mass. The coefficient Csg is 

Csg = {const.) X m,{f{^) - f{^) 

-bi^J + H^J}- ^ (17) 

As the mass difference of these two Higgs scalar masses becomes smaller, the neutron 
EDM is considerably reduced since the second Higgs scalar exchange contributes in the 
opposite sign to the hghtest Higgs scalar one as shown in the above equation. Thus, 
it is found that the second lightest Higgs scalar also significantly contributes to CP 
violation. 

In the case of tan/3 <^ 1, the contributions of Oug and O^g become very large due 
to the large lmZ2. However, these contribute to the neutron EDM in opposite signs, 
so they almost cancel each other. The remaining contribution is the one of Odg + Ogg- 

In the case of tan/3 ~ 1, the dominant contribution is the one of Odg + Ogg. In both 
regions of the large and small mH2/'mm-, the predicted neutron EDM is reduced. At 
mH2/i^m — 1) the cancellation mechanism by the second lightest Higgs scalar operates 
well, while around mH2/'iTiH3 — 0, the large mass difference of the two Higgs scalars 
leads to the small mixing between the scalar and pseudscalar Higgs bosons. 

5 Summary 

We have studied the effects of the four operators Oug, Odg + Ogg and Osg on the neutron 
EDM. The contribution of Osg dominates over that of other operators. Moreover, the 
contributions of Oug and O^g cancel out each other due to their opposite signs. This 
qualitative situation does not depend on the detail of the strong interaction hadronic 
model. Thus, the Weinberg's three gluon operator is not a main source of the neutron 
EDM in THDM. The CP violation mainly follows from the two light neutral Higgs 
scalar exchanges. Since these two exchange contributions are of opposite signs, the 



CP violation is considerably reduced if the mass difference of the two Higgs scalars is 
small. Since our predicted neutron EDM lies around the present experimental bound, 
its experimental improvement reveal the new physics beyond SM. The present upper 
hmit for is 8 x lO'^^e- cm which was given at the 26th ICHEP. Historically to reduce 
one order of magnitude for upper limit experimentally, it has taken almost 10 years. 
We hope that the rapid experimental reduction of upper limit will be performed and 
that the finite value will be reported at the close ICHEP. 
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